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■ 1 Abstract 

lO ' Most powerful jets are observed to emanate from accreting black hole systems. Recent 

quasi- stationary radiative MHD calculations reveal that jet-launching start in the 
innermost region of a transition layer (-TL), which is located between the disk and the 
^ ! corona. The plasma in the TL is dominated by virial-hot protons, advective, highly 

Q I magnetic-diffusive and rotates super-Keplerian. The amplified toroidal magnetic field 

^ ' (TMF) in the TL reach values beyond equipartition with respect to the electron 

Q ■ thermal energy, thereby considerably enhancing the Poynting energy flux and yielding 

. gravitationally unbound electron-proton outflows with a positive Bernoulli number. 

^ ! We speculate that beaming jets in GRBs may possess similar properties, pro- 

o 

^ 2 Introduction 



vided that the GRB-progenitors have rotational energy sufficient for forming a nearly 
Keplerian and massive disk. 



^ I Several observational signatures have been discovered in the last years suggesting 

^ [ that the central nucleus might be directly involved in the initiation and acceleration 

of the powerful jets observed in AGNs and BHXBs [12]. In particular, Kerr black 
holes surrounded by magnetized accretion disks are considered to be central engines 
that power the observed superluminal jets (e.g., M87 and GRS 1915+105, see jTTj . 
jHj). Those spinning black holes accreting at extremely low rates are expected to 
power electron-positron jets, whilst electron-proton jets emerge from BHs accreting 
at higher accretion rates. In the case of accretion at super-Eddington rates, the 
associated energies is likely to be trapped in the accreted matter, and subsequently 
affect the mass-growth of the central engine directly. The later case is less studied, 
and therefore it is not clear what and how the back reaction of the BH would be. 

On the other hand, recent observations reveals that GRBs might be located in 
star-forming regions, whilst their long time duration rises the possibility that they 



1 



are violent ends of massive star evolution [TH], |IH|- Recent liydrodynamical stud- 
ies of star formation, however, reveal that the first stars in the universe must have 
been massive, slowly rotating and weakly magnetized see [E], [l]- Such population of 
massive stars, the so called population III, are considered to be unstable, and they 
undergo dynamical collapse to form massive black holes pi] , PH] • Here the iron core 
of the massive star collapses to form a black hole surrounded by an accretion disk. 
Magnetic fields and neutrinos in this model are essential components for extracting 
rotational and gravitational energies to power relativistic jets along the polar axis 

In the present paper, we address the possibility that, prior to GRB-events, the 
configuration is similar to that of a BH surrounded by a massive disk, and accreting 
at extremely large super-Eddington rates. 

3 The role of magnetic fields 

In the outermost regions of a disk, we anticipate the magnetic fields to be below 
equipartition. Let r^r be a transition radius (see Fig. 1), such that for r > rtr, 
we have the usual standard disk ^j, where the ratio of magnetic to gas pressure 
P = -PmagZ-Pgas <^ 1- r < Ttr is the rcgiou where magnetic fields (-MFs) are in 
equipartiation with thermal energy of the plasma, or even stronger. In the region 
r > rtr, the Balbus-Hawley instability operates on the dynamical time scale: 
it rapidly amplifies the MP within this time scale and forces (3 to approach, but 
remains below, unity [Zj. Thus, when the matter at the outer region has completed 
one revolution around the central BH, MFs should have reached equipartition almost 
everywhere in the disk. These strong MFs will eventually suppress the generation 
of turbulence, and so terminating angular momentum transport through any sort of 
friction. 

As the rotational energy in classical Keplerian disks exceeds the thermal energy 
by at least two orders of magnitude, the generated magnetic energy via shear (mainly 
toroidal) can easily become beyond equipartition with respect to the internal energy. 

How does accretion proceed under (3=1 conditions? 

In general, accretion via strongly magnetized disks proceeds if the MFs can extract an- 
gular momentum. The efficiency of extraction depends strongly on the MF-topology. 
For example, large scale or dipolar MF-topologies are considered to be appropriate 
for extracting rotational energy from the disk or even from the hole itself and powers 
the jets. 

We may use the equations of mass and magnetic flux conservation to estimate 
the radial variation of the MFs in the vicinity of the BH. The equation of mass 
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conservation in ID reads: ^ ^ 

dt '''' r dr ^ ^ 
where S is the surface density. In the steady state case, this is equivalent to: 

r < Sf/" > -r < >^r < >= const., (2) 

where the superscripts — , + resemble the mean value of the in- and out-flowing 
plasmas, respectively. The flux r < > can be safely neglected, as most observed 
jets are expected to carry a negligibly small fraction of the matter in the disk. In this 
case, the surface density S may adopt one of the following profiles: 

r^^/^ standard disks 

r^^/^ advective disks^ (3) 

ADAF - disks^ 

Taking into account the global conservation of the magnetic flux normal to the disk, 
i.e., $ = 27rr^Sp ^ 2iTr'^Bg = const., one can easily verify that: 

standard disks 

advective disks (4) 
ADAF - disks. 

Consequently, if the disk extends down to the last stable orbit, magnetic fields 
will increase rapidly inwards, diminishing turbulence-generation and making torsional 
Alfven waves as the unigue carrier of angular momentum in the vicinity of the BH. 

Therefore, unless the imposed boundary conditions limit the amplification of 
TMF, there is no reason, why (3 should remain smaller than unity. 

Another conclusion which can be drawn from the strong dependence of Bp on 
S, irrespective of the underlying model, is that the disk must truncate near the last 
stable orbit and become a TDAT-disk (8j, otherwise a magnetic barrier will be formed 
that terminate accretion. 

To clarify the action of magnetic braking, we write the conservation equation of 
angular momentum i{= pr"^ cos 6 Q) as follows: 

Tl T2 T3 T4 





B ^ + B ^ + ——r^pu— + ^— COS 6pu— (5) 
di 39 r"^ dr dr cos 6 39 39 

TO denotes angular momentum transport via advection, Tl and T2 are for mag- 
netic extraction, and T3 and T4 are for viscous (microscopic or hydrodynamical 

"'^These solutions correspond to disk-models in which [/ « Va ~ VSs ~ 7--1/3 g^j^fj ^ jg roughly 
Keplerian JOI- here is the sound speed. 

^See ■ In writing r~^/^ we have taken into account that the vertical length scale height scales 
linearly with radius 
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Advective-Disk SS-Disk 




Figure 1: The transition from a standard to an advection-dominated disk. The 
poloidal MF (-PMF) in the outer regions is weak and below equipartition, and hydro- 
dynamical turbulence is expected to be the dominant mechanism for angular momen- 
tum transport. This is however a transient the BH-inst ability will amplify the 
MF up to equipartition on approximately the dynamical time scale. In the innermost 
region of the disk, MFs are in equipartition with the thermal energy, hence dimin- 
ishing the generation of turbulence. In this region, transport of angular momentum 
is mediated mainly by torsional Alfven waves. At i?tr we expect the equality T1=T3 
(see Eq. 5), i.e., the tranport angular momentum via hydrodynamical turbulence is 
as efficient as the vertical transport via MFs. 



turbulent) re-distribution of the angular momentum. In classical accretion disks, we 
have the equality T0=T3, whilst all other terms are neglected, traditionally. In gen- 
eral, the term Tl can be safely neglected in the equatorial near-plane if the MFs 
have a large or dipolar topology. Although it is more reasonable to consider Tl as 
the dominant term on the RHS of Eq. 5 than T3, non of the above terms will be 
neglected in the present study. 

Combining Eq.5 with the equation corresponding to the time-evolution of the 
TMF, and assuming the flow to be nearly incompressible, we obtain the torsional 
wave equation which has the approximate form: 

^ = Vl^Br. (6) 

A denotes the two-dimensional Poisson operator in spherical geometry, and Va (= 
B-p/ y/p) is the Alfven speed. The action of these waves is to magnetically brake 
the disk through transporting angular momentum from the equatorial near-plane to 
higher latitudes. Note that Bp uniquely determines the speed of propagation, hence 
the dependence of the efficiency of angular momentum transport on the -Bp— topology. 
For a nearly Keplerian-rotating medium, the removal time scale of angular momentum 
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Figure 2: The vertical distribution of rotational velocity in the disk with and without 
MFs. The action of the equipartition MPs is to transport angular momentum from the 
equatorial near-regions to higher latitudes, where the particles experience centrifugal 
acceleration and start to move outwards. 



is: 

Trem ~ pVtH/BpBt ~ r3/2. (7) 

This implies, for example, that the rate at which angular momentum is removed at 
the last stable orbit (= -Rlso) is one thousand times faster than at r = 100-Rlso- To 
avoid the innermost region of the disk from running out of angular momentum, we 
require the disk to be dynamically stable. Equivalently, the rate at which angular 
momentum is removed at any radius must be equal to the rate at which it is advected 
inwards from the outer layers, i.e., Tadv = Trom- This implies that the radial velocity 
— f/ in the disk is of the order of the Alfven speed, which is also of the order of the 
sound speed Vs, hence the terminology advection-dominated disk. 



4 Accelerating the magnetized plasma into jets 

When an isolated Keplerian-rotating particle in the equatorial near-plane is shifted to 
higher latitudes while conserving its angular momentum, its rotation becomes super- 
Keplerian locally, and therefore it starts to move outwards. Moreover, noting that 
Trem and i^disk decrease with decreasing radius, more particles in the innermost region 
of the disk are expected to be shifted to higher latitudes than in the outer regions, 
provided that /3 ~ 1 and T; ^ Tyiriai- This is because particles prefer to move along 
MF-lines than across them. As a consequence, a transition layer between the disk 
and the hot tenuous corona starts to evolve, where the matter is enforced to rotate 
super-Keplerian, and start to accelerate outwards. This configuration places an im- 
portant role for the central engine as a source for powering jets. Based on a previous 
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Radius (in 3xR5(.h units) 

Figure 3: The distribution of the velocity field in the transition region superposed 
on equally-spaced isolines of the angular velocity Q. The strong-decrease of Q with 
radius in the equatorial region relative to its slow-decrease in the TL is obvious. 

calculations it was found that that 1) the angular velocity in the transition layer 
(TL) adopts approximately the super-Keplerian profile ~ ,--5/4 2) Energy dissi- 
pation is injected primarily into the ions or protons that cool predominantly through 
fast outfiows. 3) The amplification of the generated toroidal magnetic field is termi- 
nated by the magnetic diffusion (reconnection) and fast outfiows. The width of the 
TL is pre-dominantly determined through the transverse variation of the ion-pressure 
across the jet, i.e., i^w = -Pi/VPi ~ 0.2 r, and so strongly dependent on whether the 
fiow is a one- or two-temperature plasma. In the steady-state case, this implies: 

p ~ r-^/^ Ti ~ r"^/^ U, ~ r~^/^ Bp/Bt ~ const. (8) 

We note that t/r adopts a profile and attains values similar to those in the in- 
nermost part of the disk. Provided that energy exchange between the matter in the 
disk and in the TL is efficient, the incoming matter can easily be re-directed into 
outwards-oriented motions. This implies that the Bernoulli number (Be) can change 
sign in dissipative fiows. As Fig. 4 shows. Be is everywhere negative save the TL, 
where it attains large positive values, so that the ion-plasma can start its kpc-journey. 
We note that the outfiow is sufficiently strong to shift the PMF lines outwards, whilst 
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the large magnetic diffusivity prevents the formation of large electric currents along 
the equator. In the corona however, MFs are too weak to halt the diffusive plasma 
in the dynamically unstable corona against gravity, and instead, they drift with the 
infalling gas inwards. In the case of very weak MFs (/3 < 0.1), our calculations in- 
dicate a considerably weak outflow. This is a consequence of the tendency of the 
MFs to establish a monopole like-topology, i.e, a one- dimensional MF-topology in 
which Bg — > 0. In this case, the magnetic tension becomes inefficient in feeding the 
matter in the TL with the angular momentum required for launching jets, indicating 
herewith that cold accretion disks alone are in-appropriate for initiating winds (see 

my 

Comparing the flux of matter in the wind to that in the disk, it has been found that 
Mw/ -M-d = const. ~ 1/20. The angular momentum flux associated with the wind 
is Jw/>jd = a(7Ww/ ■Mfijr^^^, where a is a constant of order unity. Consequently, 
at 500 gravitational radii almost 25% of the total accreted angular momentum in the 
disk re-appears in the wind. 

To clarify why the TL is geometrically thin, we note that in stratified dissipative 
flows the density scale height is much smaller than the scale height of the angular ve- 
locity (see P). Since the flow in TL is highly dissipative and rotates super-Keplerian, 
the time scale for torsional Alfven waves to cross the layer vertically is longer than 
the dynamical time scale at which angular momentum is advected outwards. 

On the other hand, unless there is a significant energy flux that heats up the 
plasma from below, as in the case of stars, heat conduction will always force the 
BH-coronae to collapse dynamically. To elaborate this point, let us compare the 
conduction time scale with the dynamical time scale along Bp-field at the last stable 
orbit of a SMBH: 

^ = ^ = 4.78 X 10-V,or,7,f -Ms, 

where pio, Ti lo and Aig are respectively in 10^^° g cm^^, 10^° K and in 10^ M.q units. 
This is much less than unity for most reasonable values of density and temperature 
typical for AGN-environments. In writing Eq. 4 we have taken optimistically the 
upper limit c/a/3 for the velocity, and set kq = 3.2 x 10~® for the ion-conduction 
coefficient. When modifying the conduction operator to respect causality, we obtain 
Tcond/'^dyn < Uj-/c, which is again smaller than unity. 

This agrees with our numerical calculations which rule out the possibility of outflows 
along the rotation axis, and in particular not from the highly unstable polar region 
of the BH, as ADAF-solutions predict. 
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Bernulli numbers 




Radius (in SxPg^,,., units) 



Figure 4: A snap-shot of the distribution of the Bernoulh number in two-dimensional 
quasi-stationary calculations. The decrease from large to low positive values is rep- 
resented via yellow, green and red colors. The blue color corresponds to negative 
values. The figure shows two ejected gravitationally unbound blobs of large positive 
energies in the TL. Ejection starts from the innermost region of the TL and evolves 
non- linearly. 

5 Formation of electron-proton jets: summary 

There are at least three ingredients that appears to be necessary for initiating jets 
from black hole accreting systems: 

1. large scale equipartition magnetic fields that are capable of extracting angular 
momentum from the disk efficiently, and deposit it in the viscous TL on the 
dynamical time scale, 

2. an underlying advection-dominated accretion disk for stably supplying the plasma 
in the TL with angular momentum, 

3. strong shear that is capable of generating beyond-equipartition with respect 
to the thermal energy of the electrons on the dynamical time scale. The runaway 
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Figure 5: The SED of a BH-disk-jet model of the elUptical galaxy M87. In this model, 
an accretion rate of 5 x 10~^ Eddington and T = 5 x lO^K is set to enter the domain of 
calculations through the outer boundary placed at 150 Schwarzschild radii from the 
central SMBH. The vertical scale hight of the disk at the outer boundary is taken to 
he H 0.1i?out- Additionally, a hot tenuous corona is set to sandwich the optically 
thin disk. The calculated profiles (solid lines) are superposed on the observational 
data (asterisks). The line 07 corresponds to a model in which the PMF is set to be 
in equipartition with the thermal energy, whilst TMF=0. The line 08 is similar to 
07, but the TMF is allowed to develop and reach values beyond equipartition with 
respect to the thermal energy of the electron. 

amplification of is terminated then by rapid reconnection in the magnetic- 
diffusive TL, thereby heating the proton-plasma up to the virial temperature. 

This scenario appears to apply at least for the jet in the elliptical galaxy M87 (see 
Biretta et al. 2002), where the central engine is believed to be a 3 x 10^ MIq supermas- 
sive black hole (Fig. 5). We note that since the flow in the TL is highly dissipative, we 
may expect the flow in the TL to be causally connected with the central nucleus. This 
coupling may further strengthen the role of the central BH in powering jet, therefore 
pronouncing possible observational differences between Kerr or Schwarzschild BHs.. 

6 Formating beaming jets in GRBs 

Whether GRB-progenitors are dynamically unstable massive stars or mergers of black 
holes with low mass compact stars (e.g., white dwarfs or neutron stars), the end- 
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Figure 6: The matter in the disk is advection-dominated and flow inwards with ap- 
proximately the sound speed. Only in the innermost region, where a ion-dominated 
torus is formed, matter start to move to higher latitude, while the associated angular 
momentum force the matter to start centrifugal acceleration outwards. This config- 
uration places a key role for the central engine in powering the jet in the vicinity of 
the event horizon. A spinning BH will certainly feed the jet with extra power, whilst 
a non-rotating hole will extract energy from the jet. 



product is, most likely, a BH surrounded by a massive accretion disk [?], ,2j. Unlike 
in AGNs and BHXBs, the accretion here proceeds on the dynamical time scale and 
largely exceed the corresponding Eddington rates. Therefore, almost all the energies 
associated will be trapped in the accreted matter and disappear in the hole. The 
shear in the disk excite torsional Alfven waves that carry angular momentum to 
higher latitudes, and whose speed of propagation is Va = -Bp/y/p (see Eq. 6). The 
shear in the disk amplifies -Bt in the first place, and subsequently Bp via the BH- 
instability. During the amplification time Tampi. (= Hd/VA), both Bp and B^ depend 
weakly on the density. This implies that 

^rcm - T^ampl. OC pOC M. (9) 

Consequently, the time scale for forming the TL where the jet-plasma rotate super- 
Keplerian, increases with increasing the accretion rate. 

On the other hand, since the time scale of angular momentum removal near the 
last stable orbit is of the same order as the dynamical time scale, a non-negligible 
amount of angular momentum will find its way still into the transition layer, thereby 
forming a centrifugal barrier. The accumulated matter in this region, which becomes 
extremely dense and optically thick even to neutrino emission, heats up via adiabatic 
compression, shocks, and magnetic reconnection and forms a fermi torus in which 
the effective temperatures of the electrons, protons, photons and neutrinos are of the 
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same order of the virial temperature^, i.e., 

-^e — Tp — Tj-ad — — T^virial- (-^0) 

The main coohng process here is though the loss of highly condensed matter that is 
being swallowed from below by the rapidly growing apparent horizon of the BH (see 
Fig. 6). In the case of a coUapsar, the accretion rate prior to the burst might be as 
large as 10^"^ A^Eddington Dl- 

Once the growth rate of the apparent horizon becomes comparable to, or longer 
than the centrifugal acceleration rate, the accumulated super-Keplerian fermi layer 
starts its runaway beamed explosion. 
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Discussion 



Shibata: How did you treat the corona? If the coronal density is low enough, the 
coronal MHD wind/jet would result. 



Hujeirat: Yes. The corona is hot, tenuous and without rotational support. I would 
like to stress here that obtaining coronal MHD wind, requires strong magnetic fields 
with an appropriate topology, and diffusive-free plasma. In our case, however, the 
plasma in TL is highly diffusive, inducing a rapid reconnection of the TMF, and 
heating thereby the plasma up to the virial temperature. It is because of the high 
magnetic diffusivity in the TL, torsional Alfven waves are not able to transfer angular 
momentum into the dynamically-unstable corona on the dynamical time scale. 

Fendt: You said that you obtain a jet from your calculations. But they do not look 
like that as they only show an uncoUimated wind flow along the disk surface. Please 
comment on that! 



Hujeirat: Jets are not jet-born, and coUimation is a process that occurs on later 
times and on different length scales. What interests us here is not the jet-propagation 
and collimation, but the physical mechamisms underlying their formion and accelera- 
tion. What I have shown in this talk is that in the very early phases, jet-plasmas are 
outwards-accelerated via centrifugal forces. The TL appears to be the region where 
optimal acceleration is reached: torsional Alfven waves need to propagate a small 
length scale in the vertical direction to deposit the angular momentum in the TL and 
accelerate the virial-hot protons outwards. This can be safely done on the dynamical 
time scale. Having the accelerated low-cooling protons crossed a critical radius, be- 
yond which the magnetic diffusivity becomes negligibly small, the TMF will take care 
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then to re-direct and coUimate the flow into jets. It should be noted here that the 
ion-dominated and super-Keplerian rotating plasma in the TL has a large positive 
energy, hence gravitationally unbound and potentially could expand to inflnity 
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